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Abstract 

The Debye-Hiickel theory describes rigorously the thermal equi- 
librium of classical Coulomb fluids in the high-temperature /3 — > 
regime (/? denotes the inverse temperature). It is generally believed 
that the Debye-Hiickel theory and the systematic high-temperature 
expansion provide an adequate description also in the region of small 
strictly positive values of /3 > 0. This hypothesis is tested in the 
present paper on a two-dimensional Coulomb gas of pointlike +/— 
unit charges interacting via a logarithmic potential which is equiva- 
lent to an integrable sine-Gordon field model. In particular, we apply 
a form factor method to obtain the exact asymptotic large-distance 
behavior of particle correlation functions, considered in the charge 
and number density combinations. We first determine the general 
forms of the leading and subleading asymptotic terms at strictly pos- 
itive (3 > and then evaluate their high-temperature /5 — > forms. 
In the case of the charge correlation function, the leading asymptotic 
term at a strictly positive /3 > is also the leading one in the high- 
temperature /3 — > regime. On the contrary, the /? — > behavior 
of the number density correlation function is accompanied by an in- 
terference between the first two asymptotic terms. Consequently, the 
large-distance behavior of this function exhibits a discontinuity when 
going from strictly positive values of /? > to the Debye-Hiickel limit 
/? — > 0. This is the crucial conclusion of the paper: the large-distance 
asymptotics and the high-temperature limit do not commute for the 
density correlation function of the two-dimensional Coulomb gas. 
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1 Introduction 



The object of study of the present paper is the equihbrium statistical me- 
chanics of infinite (bulk) classical (i.e., non-quantum) Coulomb fluids. For 
the sake of simplicity, we shall restrict ourselves to the case of a symmetric 
Coulomb K neutral system of two species of particles {j} of opposite 

unit charges {qj = ±1}, living in a i/-dimensional space of points r G W^. 
The system is immersed in a homogeneous medium of dielectric constant 
e = 1. The interaction energy of a set of particles {qj, rj} is given by 

E{{qj, rj}) = qjqkvilrj - rfe|), (1.1) 

j<k 

where the Coulomb potential v is the solution of the Poisson equation 

Avir) = -sj{r) (1.2) 

with Si, being the surface area of the i/-dimensional unit sphere. In particular, 

r -ln(|r|/ro) u = 2, 
v{r) = (1.3) 
[ l/|r| z/ = 3. 

The free length scale ro, which fixes the zero point of the 2D Coulomb poten- 
tial, will be set to unity for simplicity. The Fourier transform of the Coulomb 
potential given by the //-dimensional Poisson equation fll.2p exhibits the form 
l/|kp with the characteristic singularity at k ^ 0. This maintains many 
generic properties, like screening and the related sum rules [1], of "real" 
three-dimensional (3D) Coulomb fluids. For the case of pointlike particles, 
the singularity of v{r) at the origin r = prevents the thermodynamic stabil- 
ity against the collapse of positive-negative pairs of charges; in 2D for small 
enough temperatures, in 3D for any temperature. In such case the Coulomb 
potential must be regularized at short distance, e.g., by a hard-core potential 
around each particle. 

A complete analysis of the bulk Coulomb gas can be done within the 
framework of the mean-field theory of Debye and Hiickel (DH) [2j, sometimes 
called the linear Poisson-Boltzmann (PB) theory. The Debye-Hiickel theory 
describes rigorously the distribution functions of the internal Coulomb-gas 
charges in the high-temperature regime /3 — (/3 is the dimensionless in- 
verse temperature) [3]. The nonlinear ("classical" in field theory) PB theory, 
which arises as the zeroth-order term in a loop expansion of the grand par- 
tition function [1], describes rigorously the density profiles of the Coulomb- 
gas charges induced by an external charge distribution, in a specific scaling 
regime of the infinite-temperature limit [5]. The relation between the lin- 
ear and nonlinear PB theories is described in field theoretical books; see. 
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e.g., Ref. [H]. There exist many phenomenological approximations for finite 
temperatures based on heuristic extensions of the mean-field theories [7j. 

In the rigorous mathematical sense, the high-temperature — *■ regime 
is not equivalent to the case of the inverse temperature (3 being a sufficiently 
small, but strictly positive number, say [3 = 10~^^. On the other hand, it is 
generally believed that the /5 — > Debye-Hiickel theory and its improvement 
by a systematic (3 expansion adequately describe also the region of small 
positive values of /? > 0. This, at first sight natural, assumption might not 
be true for specific statistical quantities of Coulomb fiuids and therefore its 
validity has to be verified on exact results at strictly positive /3. 

To solve exactly a 3D Coulomb fiuid at strictly positive /3 is a hopeless 
task. The situation is more optimistic in the case of 2D logarithmic Coulomb 
fluids. The stability of the 2D Coulomb gas of pointlike ±1 charges against 
the collapse, associated with the 2D spatial integrability of the Boltzmann 
factor of the positive- negative pair of charges exp[/5f (r)] = |r|~^ at short dis- 
tances, is restricted to inverse temperatures (3 < 2. In this stability region the 
bulk thermodynamic properties (free energy, internal energy, specific heat, 
etc.) of the 2D Coulomb gas have been obtained exactly based on its equiv- 
alence with the (H-l)-dimensional sine-Gordon theory [8j. Later on, the 
form-factor method was applied to calculate the large- distance asymptotic 
behavior of the charge [9] and number density [10] pair correlation func- 
tions. Within the half space geometry, the surface thermodynamics (surface 
tension) of the stable 2D Coulomb gas in contact with a grounded ideal 
conductor wall was obtained through its mapping onto the boundary sine- 
Gordon model with an integrable Dirichlet boundary condition [llj. The 
exact nonperturbative asymptotic forms of the charge and number density 
profiles of the Coulomb-gas species at large distances from the conductor wall 
were calculated in Ref. [12]. It was shown that the DH theory adequately 
describes the charge profile at a small strictly positive /3, but this is no longer 
true for the number density profile. This surprising result was the primary 
motivation for the present work. 

In this paper, we reconsider the bulk 2D Coulomb gas and perform a 
more detailed form-factor analysis of the asymptotic large-distance behav- 
ior of the charge and number density pair correlation functions. In Refs. 
[HI |T0] , we have been concerned only with the exact leading large-distance 
asymptotics of the correlation functions. This information is not sufficient 
in view of present aims for which we need also terms which are subleading 
in the large-distance limit. We first review briefiy technicalities and findings 
of the previous works and then develop the form-factor formalism on higher 
levels in order to obtain the explicit forms of these subleading terms. The 
program is nontrivial and requires the application of some specific techniques 
known in the theory of integrable field models, e.g., a bootstrap procedure 
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when calculating the form-factors of higher (heavier) breathers in the parti- 
cle spectrum of the sine-Gordon theory. From the technical point of view, 
the present paper goes far beyond the previous ones [9], [10] . 

For both the charge and number density correlation functions, we first 
determine the general forms of the leading and subleading asymptotic terms 
at strictly positive (3 > and then evaluate their (3 forms. In the 
case of the charge correlation function, the leading term at positive (3 > 
is also the leading one in the /3 — >■ regime. On the contrary, the (3 
behavior of the number density correlation function is accompanied by 
an interference between the first two asymptotic terms. The large- distance 
behavior of this function therefore exhibits a discontinuity when going from 
strictly positive values of /3 > to the Debye-Hiickel limit. Equivalently, the 
large- distance asymptotics of this function at fixed temperature /3 > does 
not coincide with that obtained when the high-temperature /5-expansion has 
been performed first. This is the crucial conclusion of the paper: the large- 
distance asymptotics and the high-temperature limit do not commute for the 
density correlation function of the 2D Coulomb gas. 

The paper is organized as follows. In section 2, we introduce the notation 
and briefly summarize previous results obtained from the mapping of the in- 
finite 2D Coulomb gas onto the bulk (l+l)-dimensional sine-Gordon theory. 
The detailed form-factor analysis of the asymptotic large-distance behavior 
of the charge and number density correlation functions at a small positive 
/? > is presented in section 3. The high-temperature (3^0 behavior of 
the obtained results is discussed in section 4. A comparison is made with 
the systematic high-temperature expansion, summarized in the Appendix. 
A recapitulation and some concluding remarks are given in section 5. 



2 Sine-Gordon representation 

The bulk 2D Coulomb gas of pointlike ±1 charges is treated in the grand 
canonical ensemble characterized by the inverse temperature f3 and by the 
couple of particle fugacities z+ = z_ = z (at some places, in order to dis- 
tinguish between the + and — charges, we shall keep the notation z±). The 
grand partition function is defined by 

oo N+ N- 

E ^^e(A'..A'-). (2.1) 

where 

N 

Q{N+, N_)= Jl d'r, exp [-(3E{{qv r.})] (2-2) 
•^^ i=i 

is the configuration integral of positive and negative charges, = 
A^+ + A^_ and the interaction energy E is defined by (11.11) with v being 
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the logarithmic Coulomb potential. The infinite system is homogeneous and 
translationally invariant. Denoting the thermal average as (■ ■ ■)/3, the number 
density of particles of one charge sign g(= ±1) is defined by 



^. = (E^.,./(r-r,)) . (2.3) 
\ i I p 

Due to the charge symmetry, n+ = n_ = n/2 where n is the total number 
density of particles. At the two-particle level, one introduces the two-body 
density 

^..'(|r-r'|) = (5:5,,,/(r-r,)V,,.5(r'-rfc)\ , (2.4) 

which possesses the obvious symmetry n++ = n and n+_ = n_+. We 

consider also the pair correlation function 

V-(lr-r-l)^ -''I> -1. (2.6) 

in the charge (subscript p) and number density (subscript n) combinations: 

1 1 

g,g'=±l ^ g,g'=±l 

The 2D Coulomb gas is mappable onto the sine-Gordon theory [13]. Using 
the fact that according to Eq. (II. 2p — A/(27r) is the inverse operator of 
the Coulomb potential v and renormalizing the particle fugacity z by the 
(divergent) self-energy term exp[/?v(0)/2], the grand partition function can 
be turned via the Hubbard-Stratonovich transformation into 

„ _ JTXp exp[-^(2:)] 

^^'^ /P0exp[-5(O)]' ^'-'^ 



where 



S(z) = / d'r 



^ ~2zcos{b(f)) 



I (2.8) 



16n 

is the 2D Euclidean action of the (l+l)-dimensional sine-Gordon theory. 
Here, 0(r) is a real scalar field and / denotes the functional integration 
over this field. The one- and two-body densities of the Coulomb gas are 
expressible as averages over the sine- Gordon action (12. 8p as follows 

= ^,(e'«'"^), (2.9) 
nggi\r-r'\) = (e^'^^'^We^'''''^^^')) . (2.10) 
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The renormalized fugacity parameter z gets a precise meaning under the 
short-distance conformal normahzation 



gi60(r)g-i6<^(r')\ _ _ ^'j-^^-' as |r - r'| ^ 0; (2.11) 



for an explanation, see e.g. Ref. [1] 

The (l+l)-dimensional sine-Gordon model is an integrable field theory 
|15j . Its particle spectrum is the following. The basic particles are the soliton 
S and the antisoliton S which form a particle-antiparticle pair of equal masses 
M; for topological reason, the soliton and the antisoliton coexist in pairs. 
The S-S pair can create bound states, the so-called "breathers" {Bj;j = 
1,2, .. . < p~^}. Their number depends on the inverse of the temperature 
parameter 

The mass of the Bj breather is given by 

TJij = 2Msm(—j) , (2.13) 



and this breather disappears from the particle spectrum just when rrij = 2M, 
i.e. at the point p = Note that the breathers exist only in the stability 
region of the corresponding Coulomb gas: the lightest Bi breather disappears 
just at the collapse border p = 1 {j3 = 2), the B2 breather at p = 1/2 
((3 = 4/3), the B3 breather at p = 1/3 (P = 1), etc. 

The explicit results for the ground-state characteristics of the (1+1)- 
dimensional sine-Gordon model in the stability region were derived quite 
recently due to a progress in the method of Thermodynamic Bethe Ansatz. 
The (dimensionless) specific grand potential u, defined in the infinite- volume 
limit as 

-c^ = ^lnS, (2.14) 

was found by Destri and de Vega [16] in the form 

2 

mi 

-uj = 2.15 

8 sm(7rp) 

Here, mi is the mass of the lightest Bi breather [see formula fl2.13p taken 
with j = 1]. Under the conformal normalization (12.111) . the relationship 
between the soliton mass M and the fugacity z was established in Ref. [Ij 



r(62 



7rr(l - 62) 



2T{p/2) 



(2.16) 
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where T stands for the Gamma function. The total number density n of the 
Coulomb-gas charges is given by the standard relation 



Eqs. (I2.15P and ( \2.16h thus imply the explicit density-fugacity relationship, 
and consequently the complete bulk thermodynamics, of the 2D Coulomb 
gas in the whole stability region /? < 2 p]. 

Using the thermodynamic formulae fl2.14p - fl2.17p . the mass of the Bj 
breather fl2.13p can be expressed as follows 



Tip /Tip 

— tan — 
2 V 2 



^'^'sm(f j), (2.18) 



where k = \/2Ti^n is the 2D inverse Debye length. The high-temperature 
p ^ form of rrij reads 

lim m,=jK (j = l,2,...). (2.19) 

3 Form-factor analysis 

In a 2D integrable field theory with spectrum of particles {e} of masses 
{m^}, the two-point correlation function of local operators Oa (o is a free 
parameter) can be written as an infinite convergent series over multiparticle 
intermediate states [T5] : 

(a(r)a'(r')) = E]^ E /l^^^'^^^^'---'^^)— 



Ar=0 ■ ei,...,ejv 



'Fa,{Q^, . . . , ^i) exp ^-|r - r'| E"^., cosh^.J .(3.1) 

Here, the first = term of the series corresponds to the decoupling 
{Oa){Oa')- The next > 1 terms contain the form factors 

F,(^i,...,e^),,..,^ = (0|a(O)|Ze,(^i),...,^.^(M), (3.2) 
^--iF,,(^^,...,^i) = (Z,^(^^),...,Z,,(^i)|a'(O)|0), (3.3) 

which are the matrix elements of the local operator at the origin 0, be- 
tween an iV-particle superposition of free one-particle states, denoted as 
\Z^{9)), and the vacuum |0). The form factors, which depend on parti- 
cle rapidities {Oj} only through their differences 6jk = Oj — 9k, can be ob- 
tained explicitly in an axiomatic way [T8]. The form-factor representation of 
two-point correlation functions (13. ip is particularly useful for large distances 
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|r — r'l since it provides a systematic large-distance expansion. In the limit 
|r— r'l — * oo, the dominant contribution to the truncated correlation function 
(a(r)a'(r'))T = {Oa{r)Oa'{v')) - {Oa){Oa') comes in from a multi- 
particle intermediate state with the minimum value of the total particle mass 
SjLi^^ej! at the point of vanishing rapidities {9j — > 0}. The corresponding 
exponential decay exp(— |r — r'\J2f=i'm'tj) is modified by a slower (inverse 
power law) decaying function which particular form depends on the form 
factors. The subleading large-distance contribution to {Oa{r)Oa'{r'))T is de- 
termined by a multiparticle state with the first "excited" value of the total 
particle mass J2f=i , and so on. 

In view of the sine-Gordon representations (12.91) and (12.101) . the pair 
correlation function (12.51) is expressible as 

/ iqfe<^{r) Vfe0(r')\ 

V(k-r1)= ^ (3.4) 

This means that the local operator of interest is Oa(r) = exp[ia0(r)] with 
a = qb {q = ±1). 

The multiparticle intermediate states of the 2D sine-Gordon model can 
be created from an arbitrary number of soliton-antisoliton S — S pairs and 
breathers {Bj]j = 1,2, ... < p~^}. From the point of view of the expansion 
(13.11) . the most relevant states are those with smaller values of the total par- 
ticle mass. The general analysis of the total masses of multiparticle states 
for an arbitrary value of the temperature parameter p < 1 is complicated. 
In this paper, we are interested in the high-temperature region of relatively 
small values of the parameter p. In this region, the breathers have their 
masses (12.131) much smaller than 2M and therefore constitute the build- 
ing elements of the most relevant multiparticle states. In other words, the 
soliton-antisoliton pair of the mass 2M can be "excluded" from the construc- 
tion of the relevant multiparticle states, it only restricts the validity of the 
obtained results to specific regions of sufficiently small values of p. 

In the high-temperature regime p ^ 0, the breather masses become inte- 
ger multiples of the inverse Debye length, see formula (I2.19p . Consequently, 
a coalescence of the total masses of different multi-breather states can appear 
in that regime. In order to formalize this important phenomenon, the multi- 
breather states are grouped into subsets {Sj}'^^ such that a given subset Sj 
contains all states possessing in the regime p ^ the same total particle 
mass equal to jn. We present first few of these subsets. 

• The first subset 5*1 consists of only the one-particle Bi state with the 
lightest mass mi. The inequality mi < 2M holds in the whole stability 
interval p <1 (/? < 2). 

• The second subset 5*2 consists of the one-particle B2 state of the mass 
7712 and the two-particle -Bi-Bi state of the mass 2mi. The inequalities 
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1712 < 2mi < 2M hold in the interval p < 1/3 {(3 < 1). The equality 
1712 = 2mi takes place in the limit p 0. 

• The third subset 5*3 consists of the one-particle state of the mass m^, 
the two-particle -B1-B2 state of the mass mi +1122 and the three-particle 
BiBiBi state of the mass 3mi. The inequalities < mi + m2 < 
3mi < 2M hold for p < (2/7r) arcsin(l/3) ~ 0.216. The equalities 
m^ = mi + m2 = 3mi take place in the limit p — > 0. 

From the point of view of the next analysis, only the first two subsets are 
relevant. They contain the one-particle Bi state, the one-particle B2 state 
and the two-particle -Bi-Bi state. The exphcit form factors for these states 
are summarized in the next two paragraphs. 

For the Bi breather, the multi-breather form factors Fa{6i, . . . ,6i\f)i __i 
and ^'"^Fa'{ON, . . . , 6*1) = Fa'^O^, . . . , 6'i)i...i are presented for the exponential 
operator Oa{r) = exp[ia0(r)] in ref. [19j. Within the notation (13.21) . in the 
special case of interest a = qb {q = ±1) they read 

{0\e''''"^\Bi{e)) = -ig(7rA)^/2(e'^'^), (3.5) 
{0\e"^'"f'\Bi{e2),Bi{ei)) = -{TrX)R{ei-e2){e"^''^), (3.6) 

etc. Here, the parameter A is defined by 

A = - sin(p7r) cos (^) exp ( - / — J (3-7) 
TT V / / \ JO IT smt J 

and the function R{6) is given on the interval —271 + pir < Im(6') < — pvr by 

dt sinh(t) sinh(pt) sinh[(l + p)t] 



R{e) = A^expjs^' 



t sinh2(2t) 



X sinh 



(3. 



_ / dt sinh(t) sinh(pt) sinh((l + p)t) \ 

The function R{6) satisfies two useful relations: 

RmiO ± ivr) = ■ , . . (3.10) 

and 

Ri-e) = Siiie)Rie), (3.11) 

where 

sinh(g) + isin(p7r) 
smh(6') — 1 sm{pTT) 
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is the BiBi scattering matrix [T3]. These relations enable one to extend the 
definition fl3.8p - fl3.9l) of R{9) to arbitrary values of lm{9). For the case of 
special interest Im(^^) = 0, one gets the representation 

sinh(6') 

^ = • wm^- • f , (3-13) 



X exp < — A J 



°° dt sinh(t) sinh(pt) sinh[(l + p)t\ / 26t\ 
! sinh^(2t) V~ ; 



valid in the whole stability region < p < 1. 

The i?2 breather is a boundstate of the two Bi breathers since the BiBi 
scattering matrix f l3.12p has the B2 pole at 6' = ipn. Consequently, the one- 
particle B2 form factor can be calculated from the two-particle BiBi form 
factor fl3.6p using a bootstrap procedure |IHl [20] : 

T{0\e"''"t'\B2ie)) = ires,=o {-(7rA)i?(-ij97r + e)(e'5'<^)} , (3.14) 

where the parameter F is related to the residue of the B2 pole in the BiBi 
scattering as follows 



F = [-irese=ip.5n(^^)]'/' = ^2tan{pn). (3.15) 
Using the relation (I3.10p for 6 = —ipn + e (e ^ 0), one finds that 



\e"^''^\B2{e)) = -(vrA) 



tan(p7r)V^^ 1 



R[-i7T{l+p)] 



(e^^''^). (3.16) 



Using the explicit form factors (13.50 . (13. 6p and (13.160 in the series rep- 
resentation (13. ip for the pair correlation function (13. 4p . one finally arrives 
at 



2tan(p7r) 



VW = -gg'AKo(mir) + vrA ■^^p—(Y—yjKo{m2r) 

+ ^/(mir) + o(e-"^«^), (3.17) 



2! 

roo 



where 

^o(x)= / rie— (3.18) 

is the modified Bessel function of second kind |2T] and /(mir) denotes the 
double integral 

-00 2 J —00 2 

Based on the previous mass analysis for relevant multiparticle states, the 
asymptotic expansion (13.170 as a whole applies in the region < p < 1/3 
(0</?<l). 



10 



4 Asymptotic behavior of pair correlations 



4.1 Charge correlation function 



According to the definition of tlie cliarge correlation function hp in (12.61) . 
only the first term on the right-hand side (rhs) of Eq. (13.171) . proportional 
to qq', contributes to hp. Since the modified Bessel function Kq{x) has the 
asymptotic form 

TT \l/2 



(4.1) 



the large-distance behavior of hp, valid in the whole stability region < (3 < 
2, reads 



Kir) 



-A 



(—) 

\2mir J 



1/2 



-mir 



(4.2) 



This formula tells us that the mass mi of the lightest Bi breather is the 
renormalized inverse screening length of the charge correlation function. 
The mass mi, given by Eq. (12.181) . has the small-/3 expansion of the form 



nil 



K 



K 



sin(7r/3/(4-/3)) 



1/2 



TT 



192 



TT 

384 



(4.3) 



and the parameter A, given by Eq. (13.71) . has the small-/? expansion of the 
form 



X = /3 



1 _ (1. If^ 
\ 32 ^ 384 



237r^ 



96 2304 



13^ + 0{I3' 



(4.4) 



When /? — 0, one has nii ~ n and A ~ /3, so that the formula (14. 2 p repro- 
duces correctly the DH result (]A.9|) . We conclude that the large-distance 
behavior of the charge correlation function changes continuously when going 
from strictly positive values of /? > to the regime /3 — >■ 0, in agreement 
with the general belief. 



4.2 Number density correlation function 

According to the definition of the number density correlation function /i„ in 
(12.61) . the second and third terms on the rhs of Eq. (I3.17p . which do not 
depend on the signs of the considered charges q and g', contribute to 

hn{r) = ^^'^^^]^^p-^r^^o(m2r) + ^/(m^r) + o (e^^^) . (4.5) 

This asymptotic formula holds in the region < p < 1/3 (0 < /5 < 1). 
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The exponential decay at large r of the two terms in Eq. fl4.5p is given 
by the asymptotic of KQ{m2r) ~ exp(— m2r) and /(mir) ~ exp(— 2mir). 
At strictly positive /? > 0, the inequality m2 = 2mi cos(p7r/2) < 2mi takes 
place and therefore the large-distance behavior of hn{r) is dominated by 

h (r) ~ ^X'^-^^ i (^Y\-m.r .4g^ 

^«^^^^oo^^ 2 i?2[-ivr(l+p)] Umarj ' ^ ' 

As concerns the (3 —>■ limit of this formula, considering X ^ j3, p ^ (3/4, 
R{—i7r) 1 and m2 —>■ 2k, leads to the expression 

/in(r) ~ — — - — e (4.7) 

The large-distance form of the density correlation function is derived at 
lower orders in /3 by a systematic diagrammatic expansion in Appendix, 
see formula (lA.lSp . It is seen that the leading high-temperature term is 
of order f3^, but when one is interested in the large-distance asymptotics 
Kr ^ oo, the leading term is of order (3^. This term is twice smaller than 
the obtained result (14.71) . The reason for this inconsistency consists in the 
fact that in the regime (3 —>■ the coalescence of the inverse correlation 
lengths m2 = 2mi 2k takes place in Eq. (14. 5p . As a consequence, also the 
BiBi term (A^/2!)/(mir) in (14.50 . which is subleading for strictly positive 
f3 > 0, contributes to the large-r asymptotic behavior of hn{r) when (3 —>■ 0. 
The derivation of this additional contribution and the subsequent verification 
of the consistency of the final result with the high-temperature asymptotic 
formula (lA.lSP are the subjects of the next paragraph. 

To evaluate the /? — (or, equivalently, p — > 0) behavior of the integral 
/(mir) given by (13.190 . we first use the formula (13.130 . valid for real values 
of 9, to write down 

ft) sinh'(g) _ sm\p7c) 

p^o smh [yj + sm [pTT) cosh [U) — cos"^(]97r) 

The change of variables 6t = {6i + 62)/ 2 and 9 = 9i — d2 'i'a the integral (I3.19P 
and the subsequent integration over 9t then leads to the representation 

/(mir) ^ Kl{m,T)~ r ^ u2,!f^^''^, i^o [2mir cosh(g/2)] . 
p^o j_oo 2 cosh [9] — cos^(p7r) 

(4.9) 

It is easy to show that for any function f{9) regular at = it holds 
If co.h(^)-co.(p.) ^"'^~.?^'°'- 
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The application of this relation to the integral in Eq. fl4.9l) leads to 



2 

/(mir) ~ Knimir) - —Ko(2mir). (4.11) 
p^o 2 

Thus, 

^I{m,r) = ^KliKv) - ^-^K,{2Kr) + 0{(3% (4.12) 
At asymptotically large distance, 

^I{m,r) ^e-^--^-^(-^y^\-^- + 0{P% (4.13) 
2! ^ ^ Anr 16 KAnrJ ^ ^ ^ 

Summing up the rhs of this formula with the rhs of the previous formula 
(14.71) . one recovers correctly the asymptotic result (lA.lSP obtained by the 



systematic (3 expansion. We conclude that the large-distance behavior of the 
number density correlation function undertakes an abrupt change, namely 
the discontinuity, when going from strictly positive values of /3 > to the 
/9 — >■ regime. 

We have shown by the exact calculation that the large- distance asymp- 
totics of the density correlation function at fixed temperature /? > does not 
coincide with that obtained when the high-temperature /5-expansion has been 
performed first. This is equivalent to saying that the large- distance asymp- 
totics and the high-temperature limit do not commute for this function. Such 
phenomenon is in contradiction with the "usual" physical intuition. 



5 Conclusion 

In the present paper, we took advantage of the exact solvability of the bulk 
2D Coulomb gas to study the large-distance behavior of correlation functions 
between charged particles. Using the form-factor technique for the equivalent 
(l-l-l)-dimensional sine-Gordon theory, we have expressed in Eqs. fl3.17p - 
(I3.19P the leading and subleading asymptotic terms of particle correlation 
functions in terms of the masses of breathers belonging to the sine-Gordon 
particle spectrum. 

The result for the charge correlation function hp{r) fl4.2l) - fl4.4p . valid in 
the whole stability region of the Coulomb gas < (3 < 2, has the generally 
anticipated property: in the /? — regime, it reduces continuously to the 
DH result (lA.9p . This means that heuristic extensions of mean-field theories 
to finite temperatures are reasonable when they are based on the charge- 
charge correlations. 

On the other hand, the formula fl4.6p for the asymptotic decay of the 
number density correlation function /i„(r), valid in the region < /3 < 4/3, 
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does not reproduce in the /3 — regime [see Eq. fl4.7p ] the result of the 
high-temperature expansion (lA.lSp . The reason for this inconsistency con- 
sists in the fact that when j3 —>■ the term (A^/2!)/(mir) in Eq. (14.51) . 
which is sub leading for strictly positive (3 > 0, interferes with the leading 
one and also contributes to the asymptotic result. Taking into account the 
asymptotic formula (14.131) for this subleading term, one recovers correctly 
the high-temperature formula flA.lSp . As a consequence of the above sce- 
nario, the large-distance behavior of the number density correlation function 
undertakes a discontinuity when going from strictly positive values of (3 > 
to the /3 — i> regime. The high-temperature expansion (lA.lSP therefore does 
not reflect adequately the large-distance behavior of the number density cor- 
relation at strictly positive /3 > 0. This phenomenon contradicts the general 
belief and one has to be very careful when extending the DH description of 
number density correlations to finite temperatures. 

We notice that the derivation of the standard Debye-Hiickel theory is 
based on electrical quantities and it is perhaps not surprising that the be- 
havior of quantities related to the number density is not always adequately 
reproduced in this high-temperature theory. 

The anomaly in the large-distance behavior of the number density corre- 
lation function hn{r) was observed due to the availability of the exact (and, 
therefore, nonperturbative) description of the 2D Coulomb gas. The anomaly 
could be observed perturbatively only after the resummation of specific dia- 
grammatic contributions in all /3 orders of the large-distance decay of hn{r). 
It would be interesting to reveal the resummation mechanism because the 
described anomaly might be present in 3D Coulomb fluids, too. 



Appendix: Diagrammatic expansion 

In this appendix, we derive the asymptotic large-distance form of the charge 
and number density correlation functions in the bulk 2D Coulomb gas, at 
lower orders in (3. The f3 expansions of the correlation functions must be 
taken for a fixed value of the inverse Debye length n = ^/2^^J3n, which only 
sets the length scale. 

For the considered Coulomb gas with the charge symmetry, the ordinary 
Ornstein-Zernike (OZ) equation splits into two independent relations for the 
charge and density functions [22] 

hp = Cp + Cp*n* hp, (A.l) 

hn = Cn + Cn*n*hn, (A.2) 

where * denotes a convolution product and the charge and density direct 
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correlation functions are defined in analogy with Eq. fl2.6p as follows 

^ q,q'=±l q,q'=±l 

In the renormalized Mayer expansion of the excess Helmholtz free energy 
in density [22], the chains of simple —jSv bonds are summed up into the 
renormalized bonds K, defined implicitly by the relation 

K = -I3v + {-I3v) *n*K (A.4) 

with V being the Coulomb potential. In the infinite 2D space, one has 

K{r) = -iSKoinr), (A.5) 

where Kq is the modified Bessel function of second kind (13.181) . 

In the renormalized-bond format, the charge direct correlation function 
is expressible as 

00 

c,(r) = -Mr) + 5:/5^c(/-)(r), (A.6) 
i=3 

where only such renormalized graphs of the excess Helmholtz free energy 
contribute to the coefficients {c^J^}'^^. which have their two root vertices 
with an odd bond-coordination and their field vertices with an even bond- 
coordination. In particular [9], 

(r) = -Ik'oM, (A.7) 
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c(^)(r) 



^Ko{Kr) [ d\Kr')K',{Kr')K^,{K\r - r'|), (A.^ 



etc. Inserting the leading Cp(r) = — /9t>(r) into the OZ relation ( lA.ip . one 
gets for hp nothing but the definition (]A.4I) of the renormalized bond K. 
Therefore, at the lowest order in /3, 

hp{r) = -PKoinr) ~ -fi ( — ) exp(-«:r) (A.9) 

since -^'0(2;) has the asymptotic form (14. ip . It was shown in Ref. [9] that the 
consideration of the terms j = 3, 4 in the expansion (\A.6\i with the corre- 
sponding coefficients (]A.7p and (lA.SP implies the large-distance asymptotic 
behavior (14. 2 p with the f3 expansions of the parameters mi (14. 3p and A (14. 4p 
up to the indicated order. 

Within the renormalized-bond formalism, the number density direct cor- 
relation function is given by 

c^{r) = -K\r) + J2P^c^^\r). (A.IO) 



15 



Here, the leading term of order (3'^ corresponds to the renormahzed Meeron 
(watermelon) diagram and renormalized graphs of the excess Helmholtz free 
energy, contributing to the coefficients {c^^''}°f,4, are the ones which have their 
two root vertices as well as field vertices with an even bond-coordination. 
Inserting the leading c„(r) = ^K'^{r) of order into the OZ relation flA.2l) . 
the convolution term is easily seen to be of higher order Thus, at lowest 
order in /3, /i„(r) = c„(r), 

K{r) = ^KliKv) ^r^J^-eM-'^^r). (A.ll) 

Since the sum on the rhs of flA.lOl) starts from j = 4, the f3^ term of hn has 
its origin exclusively in the convolution term of the OZ relation flA.2p taken 
with c„(r) = hn{r) = K'^{r)/2\. Consequently, 

h^ir) = ^Klinr) + ^J{Kr) + 0(/5^), (A.12) 
where the integral J is defined by 

J{r) = j'^Kl{v')Kl{v-r'). (A.13) 
In terms of the 2D Fourier transform of Kq{t) 



dV , In (k/2) + Jl + ik/2) 



27r 

J is expressible as follows 



J(r) = / ^e^'^-G^k). (A.15) 
Let us put r = (0, r) in the Cartesian notation. Since the function G'^{k) has 



simple poles at ky = ±1^4 + /c^, the integration over the vector component 
ky can be performed explicitly as the contour integration in the complex 
plane by using the residuum theorem, with the result 

J{r) = — / ^ , exp (-r^/IT¥) . (A.16) 



lA + kl 

The integral in flA.16D is equal to -ft'o(2r) and so J(r) = [tt"^ /A)KQ{2r). Eq. 
()A.12p thus takes the form 

K{r) = ^KUkt) + ^-^K,{2Kr) + 0{(3'). (A.17) 
2 io 
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At asymptotically large distance, 



K{r) ~ ^e-^- + ^(-^) e-'-^ + 0{(3% (A.18) 
r^oo 4f^j- 16 \4KrJ 

Here, the leading high-temperature term is of order But when one is in- 
terested in the large-distance nr ^ oo asymptotics of the density correlation 
function, like in the present paper, the leading term is of order (3^. 
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